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Abstract 

We study an approximate version of the Schwinger-Dyson equation that controls the nonpertur- 
bative behavior of the ghost-gluon vertex, in the Landau gauge. In particular, we focus on the form 
factor that enters in the dynamical equation for the ghost dressing function, in the same gauge, 
and derive its integral equation, in the "one-loop dressed" approximation. We consider two special 
kinematic configurations, which simplify the momentum dependence of the unknown quantity; in 
particular, we study the soft gluon case, and the well-known Taylor limit. When coupled with 
the Schwinger-Dyson equation of the ghost dressing function, the contribution of this form factor 
provides considerable support to the relevant integral kernel. As a consequence, the solution of 
this coupled system of integral equations furnishes a ghost dressing function that reproduces the 
standard lattice results rather accurately, without the need to artificially increase the value of the 
gauge coupling. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 
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I. INTRODUCTION 



One of the few nonperturbative frameworks available for the study of the infrared sector 
of QCD in the continuum are the Schwinger-Dyson equations (SDEs), which govern the 
dynamics of the basic Green's functions of the theory [l-4]. Despite the well-known limita- 
tions intrinsic to this formalism, a variety of theoretical and technical advances have provided 
new valuable insights on some of the most fundamental nonperturbative phenomena of QCD, 



such as quark confinement, dynamical mass generation, and chiral symmetry breaking 



Particularly important in this ongoing effort is the systematic confrontation of the SDE 
predictions with the results of large-volume lattice simulations js-Jjj], leading not only to 
quantitative refinements, but, in some cases, to critical revisions of the underlying physical 
concepts j?! Il5- 17]. 

The quantitative understanding of the ghost sector of QCD constitutes a long-standing 
challenge for the SDE practitioners. Without a doubt, the most fundamental quantity 
in this context is the ghost propagator, D(p 2 ), and the corresponding dressing function, 
F(p 2 ) = p 2 D(p 2 ); in fact, the infrared behavior of the latter, in the Landau gau ge ( LG), has 



18 



2l|. 



been traditionally associated with a particular realization of color confinement 

In recent years, various lattice studies, both in SU (2) and SU(3), together with numerous 
analytic approaches, find a massless ghost propagator with an infrared finite dressing func- 



tion 



QflHHI 



|22J. In addition, in the same gauge, the gluon propagator obtained on the 
lattice is finite in the deep infrared, supporting the notion of an effectively massive gluon. 
In fact, the dynamical gluon mass generation, first proposed in 23], and further developed 
in a number of recent works, provides a unified explanation for the observed finiteness of 
both aforementioned quantities Specifically, an infrared finite F(p 2 ) emerges as 

a direct consequence of the massiveness of the gluon propagator: such a gluon propagator, 
when inserted in the SDE of the ghost propagator, saturates the logarithms associated with 
the F(p 2 ), thus making it finite at the origin. 



f-consistent SDE analysis is 
2811 . In fact, even when one 



However, what has been more difficult to obtain from a se 
the eath-e shape a*d s,e of ^provided* the lattice g 
substitutes into the ghost SDE the gluon propagator furnished by the lattice, but keeping 
the ghost-gluon vertex at its tree-level value, the resulting F(p 2 ) is significantly suppressed 
compared to that of the lattice to reproduce the lattice result, one has to artificially 
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increase the value of the gauge coupling from the correct value a s = 0.22 to a s = 0.29 [29 ]. 

It would seem, therefore, that the main reason for the observed discrepancy ought to 
be traced back to the way in which the fully dressed ghost-gluon vertex, T u , appearing in 
the ghost SDE, is approximated. Even though preliminary lattice studies indicate that the 
deviations of T v from its tree-level value are relatively moderate 
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34|, the highly non- 



linear nature of the ghost SDE may lead to considerable enhancements. In fact, a modest 
increase of the relevant vertex form factor in the region of momenta that provide the largest 
support to the ghost SDE may account for the bulk of the required effect. 

The purpose of this article is to obtain a reliable "first-principle" approximation for this 
important missing ingredient. Specifically, we will determine the relevant vertex form factor 
from an approximate version of the SDE satisfied by the vertex Y v itself, in the LG. To be 
sure, the vertex SDE has a complicated skeleton expansion, involving various unknown (or 
only partially known) quantities, such as multiparticle kernels. The basic approximations 
we employ at the level of the vertex SDE are: (i) we consider only the first two diagrams 
in this expansion; this corresponds to the "one- loop dressed" truncation [35], and (ii) inside 
these diagrams we replace full vertices by their tree-level values, but keep fully dressed ghost 
and gluon propagators, (iii) for the numerical analysis of the resulting integrals, we use as 
input for the full gluon propagators the lattice data of 11]. 

The tensorial decomposition of Y v consists of two form factors [see Eq. (12.21) ]; however, 
given that this vertex will be inserted in the ghost SDE, written in the LG, only the co- 
factor A(— k, — p, r) of the ghost momentum p v survives. In the present study we deter- 
mine A(—k, —p, r) for two particular kinematic configurations, soft gluon (k — > 0) and soft 
ghost (p — > 0), thus converting it, in both cases, to a function of a single momentum only, 
A(0, —p,p) and A(—k, 0, k), respectively. In fact, as we will explain in detail in Sec. IHIt the 



case where p — > is equivalent to the standard Taylor limit [36|, |37|. 
Our main results may be summarized as follows. 

(i) In the soft gluon limit, the result obtained for A(0, —p,p) displays a moderate peak 
around 1 GeV, corresponding to a 20% increase with respect to the tree-level value; this 



result compares rather well with the existing lattice data |32l. I33j]. Of course, this particular 
kinematic configuration is not relevant for the ghost SDE, but serves as a preliminary test 
of the overall faithfulness of the approximations employed. 

(ii) The numerical solutions for the coupled system of integral equations determining 
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F(p 2 ) and A(—k, 0, k) gives rise to a ghost dressing function that is in excellent agreement 



with the lattice data 



The corresponding solution for A{— k, 0, k) is characterized by a 



rather pronounced maximum, centered again around 1 GeV, reaching a value of about 1.5. 
In this analysis we use a s = 0.22, which corresponds to the momentum-subtraction (MOM) 



value for the point \x = 4.3 GeV 38|, used to renormalize the gluon propagator obtained 
from the lattice. 

The article is organized as follows. In Sec. [TT] we introduce the necessary notation, and 
set up the SDE for the ghost dressing function, paying particular attention to the way that 
the fully dressed ghost-gluon vertex enters in it. In Sec. II III we carry out the analysis at the 
level of the SDE of the ghost-gluon vertex, and derive the corresponding closed expressions 
in the two kinematic limits of interest. In Sec. IIVI we present the numerical treatment of 
the equations derived in the previous sections. In particular, we first compute the case of 
the soft gluon, and then we proceed to the solution of the coupled system. Finally, our 
conclusions and discussion are presented in Sec. |Vl 



II. GHOST DRESSING FUNCTION AND THE GHOST-GLUON VERTEX 

In this section we introduce the SDE for the ghost propagator in the LG, and discuss some 
of its basic properties and features. Of particular interest is the dependence of this equation 
on the surviving component of the ghost-gluon vertex, and the numerical implications of 
approximating it by its tree-level value. 

Our starting point is the full ghost-gluon vertex, shown in Fig. (TJ and denoted by 

r:"(-k,-p,r)= g f" t T,(-k,-p,r), r = k + p, (2.1) 



t* 



k + p 



P 



FIG. 1: The fully dressed ghost-gluon vertex. 
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r„(-fc, -p, 



k+ P ) 




FIG. 2: The SDE for the ghost propagator given by Eq. (|2.T|) . The white blobs represent the fully 
dressed gluon and ghost propagators, while the black blob denotes the dressed ghost-gluon vertex. 

with k representing the momentum of the gluon and p of the anti-ghost. The most general 
tensorial structure of this vertex is given by 

T u (-k, -p, r) = A(-k, -p, r) p v + B(-k, -p, r) k v ; (2.2) 

at tree-level, the two form factors assume the values A^(—k,—p,r) = l and 
B^(—k, —p,r) = 0, giving rise to the bare ghost-gluon vertex r{? = p v - 

The form factors A and B may be formally projected out by contracting T v with the 
vectors 

A/. \ _ k2 Pu ~ (k-p)k v _ p 2 k v - (k ■ p)p u 

£ u{k,p)- k2p2 _ {k . p)2 > e u {k,p)- k2p2 _ {k . p)2 , V-V 

namely 

A(-k, -p, r) = e A u (k, p)V(-k, -p, r) , B(-k, -p, r) = e B u (k, p)Y v {-k, -p, r) . (2.4) 

Of particular importance for the analysis that follows is the so-called "Taylor limit" of 
the ghost-gluon vertex, corresponding to the case of vanishing ghost momentum, r = 0, 
p = —k. In this special kinematic configuration, the T u (—k, —p,r) of Eq. ( 12. 2ft becomes 

T v (—k, k, 0) = -[A(-k, k, 0) - B(-k, k, 0)}k u . (2.5) 

Closely related to this limit is the well-known Taylor theorem, which states that, to all orders 
in perturbation theory, 

A(—k, k, 0) - B(—k, k, 0) = 1; (2.6) 

as a result, the fully-dressed vertex assumes the tree- level value corresponding to this par- 
ticular kinematic configuration, i.e., T u (—k,k,0) = —k u . 
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After these introductory comments, let us turn to the SDE for the ghost propagator, and 
examine in some detail how the ghost-gluon vertex affects its structure. The relevant SDE 
is diagrammatically represented in the Fig. |2J Using the momenta flow and Lorentz indices 
as indicated in Fig. [2J the ghost SDE can be written as 

l D-\p 2 )= ip 2 -g 2 C A [ r^{k,-k-p,p)A^{k)T u {-k,-p,k + p)D(k + p), (2.7) 



k 



where Ca denotes the Casimir eigenvalue of the adjoint representation (N for SU(N)), 
d = 4 — e is the space-time dimension, and we have introduced the integral measure 

with \i the 't Hooft mass. In the LG, the gluon propagator A M „(g) has the transverse form 

A^(q) = -iP^{q)A{q 2 ) , (2.9) 

with 

PM = 9^-^, (2.10) 

the usual projection operator. 

Clearly, due to the full transversality of A AtI/ (fc), any reference to the form factor B 
disappears from the ghost SDE of Eq. (12. 7p . Specifically, substituting Eq. (12. 2p into Eq. (12. 7p 
we obtain 



F- L (p 2 ) = l + tg 2 C A 



(k-p) 



A(-k,-p,k+p)A(k)D(k + p), (2.11) 



k L k2 P 2 

where we have introduced the ghost dressing function, F(q 2 ), defined as 

D(q 2 ) = ^. (2.12) 

The renormalization of Eq. (12.111) proceeds through the replacements 

A R (q 2 ) = Z^Aiq 2 ), 

F R (q 2 ) = Z: l F(q 2 ), 

r fl(<LP,0 = Z l T v (q,p,r), 

g R = Z g l g = Z^Z)[ 2 Z c g, (2.13) 

where Z&-, Z c , Z\, and Z g are the corresponding renormalization constants; the dependence 
of the above quantities on the renormalization point [i is suppressed. In the MOM scheme, 



usually employed in the SDE analysis, the renormalization conditions imposed are that at fi 
the corresponding Green's functions assume their tree-level values, e.g., A~ x (g 2 = /i 2 ) = fi 2 , 



and F R (q 2 = fi 2 ) = 1 37]]. Note also that, in the LG, the form factor A is ultraviolet finite 
at one-loop, and therefore, no infinite renormalization constant needs to be introduced at 
that order for T v . In fact, one usually invokes Taylor's theorem [see Eq. f)2.6p ]. in order to 
finally set Z\ = 1 to all orders (see discussion in Sec. IHIj) . 
Then, the SDE becomes 

(k-p) 2 ' 



F- 1 (p 2 ) = Z c + zg 2 C A [ 

Jk 



A(-k,-p,k + p)A(k)D(k + p), (2.14) 



k 2 p 2 

where we have suppressed the subscript "R" to avoid notation clutter. The actual closed 
expression of Z c is obtained from Eq. (I2.14p itself, by imposing the aforementioned MOM 
renormalization condition on F~ 1 (p 2 ). 

Evidently, the explicit dependence of Eq. (12. 14ft on A(—k, —p, k + p) requires the use 
of the corresponding vertex SDE, thus converting the problem of determining F(p 2 ) into a 
coupled SDE system. The usual way to circumvent this technical complication has been 
to simply approximate A(—k, —p, k + p) by its tree-level value, setting into Eq. (12 . 14j) 
A(—k, — p, k + p) = 1. 

Then, after proper renormalization along the lines discussed above, and passing to the 
Euclidean space following the standard rules, one solves Eq. (12.141) numerically, using the 



lattice data of 11] as input for the gluon propagator. Note that this latter propagator is 
renormalized within the MOM scheme, by imposing the standard condition A -1 ^ 2 ) = /i 2 
at fi = 4.3 GeV, namely the deepest available point in this set of lattice data; then, the 
corresponding value for a s = g 2 /47r that one should use is a s (4.3GeV) = 0.22. However, for 
this particular value of a s , the solution obtained from Eq. (12.141) lies considerably below the 
lattice data for F(p 2 ), as can be clearly seen from the (blue) dotted curve of Fig. [3j In order 
to obtain a close coincidence with the lattice, one must increase the value of a s (4.3GeV) to 
0.29, thus obtaining the (red) continuous curve in Fig. [3j 

It is, of course, natural to attribute the observed discrepancy to the aforementioned simple 
approximation employed for the ghost-gluon vertex. Therefore, to ameliorate the situation, 
we will determine this form factor from its corresponding SDE, in a certain kinematic limit 
that is relevant for the situation at hand. Specifically, given that A(—k, —p, k + p) is a 
function of three variables, p 2 , k 2 , and the angle between the two (appearing in the inner 




3.2- 



Ghost dressing function 

a Lattice V=64 4 and ()=5.7 
* Lattice V=80 J and (1=5.7 
Ghost SDE (bare vertex) for a(n)=0.22 



0.8 



1E-3 



0.01 



0.1 



10 



100 



p 2 [GeV 2 ] 



FIG. 3: Comparison of the ghost dressing function, F(p 2 ), obtained as solution of the ghost SDE 
when the ghost-gluon vertex is approximate by its bare value, with the lattice data of Ref. fll| . 
The (red) continuous curve represents the case when a s (4.3GeV) = 0.29 whereas the (blue) dotted 
curve is obtained when a s (4.3GeV) = 0.22. 

product p ■ k), a full SDE treatment is rather cumbersome, and lies beyond our present 
technical powers. Instead, we will consider the behavior of A(— k, —p, k + p) for vanishing 
p\ to that end, we start out with the Taylor expansion of A(— k, —p, k + p) around p = 0, 
and we only keep the first term, A(— k,0,k), thus converting A into a function of a single 
variable. 

We emphasize that the limit p — > is taken only inside the argument of the form factor 
A, but not in the rest of the terms appearing in the SDE of Eq. ( I2.14I) . Specifically, following 
the procedure explained in detail in the next section, one isolates from the ghost-gluon SDE 
the contribution proportional to p u , taking the limit p — > in the accompanying scalar 
co-factor, thus arriving at a form T u (— k, —p, k + p) = p u A(—k, 0, k). Equivalently, in terms 
of the projectors introduced in Eqs. ( 12. 3ft and ( 12.41) . one has 



A(-k,0,k) = lim{eC(k,p)T u (-k,-p,k + p)} . 



(2.15) 



Thus, the approximate version of the SDE in Eq. ( 12.141) reads 



F- 1 (p 2 ) = Z c + tg 2 C A [ 1- 



Jk 



(k-p) 
k 2 p 2 



A(-k,0,k)A(k)D(k + p). 



(2.16) 



S 




FIG. 4: (A) The complete SDE of the ghost-gluon vertex. Notice that we have set up it with 
respect to the anti-ghost leg. (B) Diagrams included in the skeleton expansion of the ghost-gluon 
kernel that we will consider in our analysis. 

III. THE GHOST-GLUON VERTEX 

In this section we derive in detail the nonperturbative expression for the form factor A, 
in two special kinematic configurations: (i) the soft gluon limit, in which the momentum 
carried by the gluon leg is zero (k — 0), and (ii) the soft ghost limit, where the momentum 
of the anti-ghost leg vanishes (p = 0). 

A. General considerations 

The starting point of our analysis is the SDE satisfied by the ghost-gluon vertex, whose 
diagrammatic representation is shown in panel (A) of Fig. HI One observes that the relevant 
quantity, which controls the dynamics of this SDE, is the four-point ghost-gluon kernel. 
For the ensuing analysis we will carry out the following main simplifications: 
(i) The ghost-gluon kernel will be replaced by its "one-loop dressed" approximation; 
specifically, in the corresponding skeleton expansion we will only include the diagrams ap- 
pearing in panel (B) of Fig. HI Thus, the approximate version of the SDE that we employ 
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may be cast in the form 

T v (-k, -p, k+p) =p u - -g 2 C A [(di) u - (d 2 )u\ , (3.1) 
where the diagrams (di) are given by 

{d!) v = !rf^{l)T a D{l + k+ P )T v D{l+p), 



(d 



2 v 



Jvf^{l)V vaa ^{l - k)TpD(l+p) . (3.2) 



For notational simplicity, we have suppressed the arguments of the momenta in all ver- 
tices; the latter may be easily recovered from the figures and the conventions established in 
Sec. HIl Note also that, in the LG that we use, the gluon propagators appearing in the above 
expressions assume the completely transverse form of Eq. (I2.9p . 

(ii) The (multiplicative) renormalization of Eq. (13. ip proceeds in the standard way. 
Specifically, in addition to the renormalization constants and relations given in Eq. (I2.13p . 
one must introduce the vertex renormalization for the three-gluon vertex, to be denoted by 
Z3, namely T v ° a = Z^T uaa , together with the corresponding relation for the coupling renor- 
malization, namely g R = Z^Z^ 2 g. From this relation, and the last of Eq. (12. 13ft . one has 
that Z^ 1 Za = Z^ X Z C . Then, it is straightforward to show that the contributions of g 2 {d\) v 
and g 2 {d\) v maintain the same form after renormalization; in fact, this property may be eas- 
ily established by grouping the integrands in terms of the standard renormalization-group 



invariant quantities formed by (gY^^^D) and (gT uaa A 3 ^ 2 ) 39[. Thus, the renormalized 
version of Eq. (13 .ip reads 

Tl(-k,-p,k+p) = Z 1 L v -^£C A [(di)£ -(<*»)£]} , (3.3) 

where the Z\ comes from the renormalization of the T u (—k, —p, k + p) on the lhs. 

In what follows we will set Z\ = 1. In the case of the soft ghost configuration, p = 0, 
(which, as we will see, is equivalent to the Taylor kinematics), this choice is imposed by 
Taylor's theorem, see Eq. (I2.6p . On the other hand, in the case of the soft gluon configuration, 
k — 0, this choice constitutes an approximation, in the sense that it is motivated by the 
one-loop finiteness of the (LG) T u , but is not enforced by an analogous all-order relation. 

(iii) In the two aforementioned diagrams, (di) and (d 2 ), we will keep fully dressed prop- 
agators, but will replace the fully dressed three-gluon vertex appearing in graph (rf 2 ) by the 
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corresponding tree-level expression, namely 

F a »v(q, r, p) -»■ rg, r, p) = (r - p) a ^ + (p - g) M flW* + (q - r) v g ail . (3.4) 

Furthermore, as will be explained in the corresponding subsections, additional approxima- 
tions will be imposed on the fully dressed ghost-gluon vertices, depending on the specific 
details of each kinematic case considered. 

B. Soft gluon configuration 

We begin with the analysis of the soft gluon configuration, k = 0. Evidently, in this case 
the ghost-gluon vertex becomes a function of only one momentum, p, and may be described 
in terms of a single form factor, namely, 

I\,(0, -p,p) = A(p)p u ; A{p) = A(0, -p,p) . (3.5) 

Therefore, setting k = in Eq. (I3.ip . one is able to isolate the form factor A by means of 
the projection 

A(p) = l- t -g 2 C A [(d 1 )-(d 2 )]; = 2 = 1,2, (3.6) 

L p 

where the diagrams (di) are obtained from those of Eq. (I3.2p in the limit k — > 0. 

The particular kinematic configuration considered here allows one to derive a linear in- 
tegral equation for the unknown quantity A(p). This becomes possible because, in the limit 
k = 0, the vertex T v entering in graph (d\) becomes T u (0, —l—p,l+p). Thus, the integral 
(di) contains A(0, —l—p,l + p), giving rise to an integral equation for A(0, — p, p). Unfor- 
tunately, this favorable set of circumstances does not apply to the remaining ghost-gluon 
vertices, namely, IV and in graphs (di) and (c?2), respectively; their arguments depend 
on all possible kinematic variables, and the inclusion of the full A would give rise to a (non- 
linear) integral equation, too complicated to solve. We therefore approximate all remaining 
ghost-gluon vertices by their tree-level expressions. 

After these comments, and use of the notation introduced in Eq. (13.61) . the diagram (d±) 
reads 

(di) = l yl'p^ -P) 2 ~ lV]D 2 (l)A(l+p)A(l) . (3.7) 
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To evaluate the contribution of diagram (<i 2 ) notice that, with the gluon propagators in 
the LG, and the bare three-gluon vertex of Eq. (|3.4p . we have that 



P^(/)P^(/)41(0, l, -l) = 2l v P>>P{l) . 



(3.8) 



Applying this result we get 



(4) = 2 l^iiy-(l-pf]A\l)D(l + p). 



(3.9) 



The final answer is obtained by substituting Eq. (13. 7p and Eq. (I3.9P into Eq. (13.61) ; it will 
be written directly in Euclidean space, using the standard transformation rules, 



- q 



and setting 



A E (q 2 E ) = -A(-ql)- D E (q 



-D(-q 



(3.10) 



l 2 = t; p 2 = x; (l+p) 2 = z; (I ■ p) = \/xi cos 6 ; 

dtt / d9sm 2 9 



, lE j (2vr)4 (27r)3 7 J 
Specifically (we suppress the subscript "E"), 

_. /^i poo 

a s C A 



(3.11) 



A(x) 



An 2 



° [A(z) 
dtVxiF 2 (t)A(t) /^sin 4 ^cos^ 

Jo L z 



a s C A 



2vr 2 



/•OG P7T 

/ dtVxitA 2 (t) / ^sin 4 ^cos^ 
Jo Jo 



F(z) 



(3.12) 



where we have used g 2 = 4na s , and Eq. ( 12 . 1 2 j) in order to express the ghost propagators in 
terms of their dressing functions. 

Notice that, in the limit x — 0, namely when the momentum of the ghost leg is also zero, 
we recover from Eq. (13. 12ft the tree-level value of the form factor, i.e., A(0) = 1. 



C. Soft ghost configuration (Taylor kinematics) 

We next turn to the case that, according to the discussion presented in Sec. [Til is expected 
to improve the treatment of the ghost SDE. Specifically, in this subsection we will derive an 
approximate version for A in the soft ghost configuration, to be denoted by 



lim A(-k, -p, k+p) = A(-k, 0, k) = A(k) 

p->0 



(3.13) 
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However, before proceeding to this derivation, we will demonstrate that the form factor 
A(—k, 0, k) obtained in the soft ghost configuration is none other than the form factor 
A(— k, k, 0), appearing in the constraint imposed by Taylor's theorem, given by Eq. (12. 6p . 
To prove that, let us rewrite the SDE of the ghost propagator, Eq. (12. 7p . dressing this time 
the left ghost-gluon vertex instead of the right, i.e., 

iD-\p 2 ) = tp 2 -g 2 C A [ T^k,-k-p,p)A^(k)T^{-k,-p,k+p)D{k+p), (3.14) 

Jk 

where we have maintained the same momenta flow and Lorentz indices as in Fig. [2j There- 
fore, using Eq. (Q for the T M in Eq. (15TT4"]) . we get (in the LG) 

F-\p 2 ) = l + ig 2 C A [ 

Jk 

Evidently Eqs. (12. lip and ( I3.15P must furnish an identical result for F(p 2 ), since the answer 
cannot depend on which of the two vertices one chooses to dress. Thus, the form factor A 
is forced to satisfy the equality 

A(-k, -p,k + p) = A{k, -k-p,p). (3.16) 

Given that, due to Lorentz invariance, the dependence on the momenta is quadratic, i.e., 
A(k 2 ,p 2 , k 2 + p 2 + 2k ■ p), we have immediately that 

A(k, -k-p,p) = A(-k,k+p, -p). (3.17) 

So, combining Eqs. (I3.16P and (13.171) . we arrive at the relation 

A(-k, -p, k + p) = A(-k, k + p, -p), (3.18) 

which states that, in the LG, the form factor A of the gluon-ghost vertex is invariant under 
the exchange of the momenta of the ghost and anti-ghost legs. Notice that this invariance is 
known to be a consequence of a global SL(2, R) symmetry between the ghost and anti-ghost 
fields which implies that the LG is a ghost-anti-ghost symmetric gauge fixing choice. 
Finally, setting p = in Eq. (I3.18p . one obtains the announced result, that is, the A obtained 
in the soft ghost limit coincides with that of the Taylor kinematics. 

As mentioned above, the fact that the kinematic situation considered here is equivalent 
to the Taylor limit, imposes, in a natural way, the value Z\ = 1 for the renormalization 
constant appearing in Eq. (13.31) . 
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(k-pY 
k 2 p 2 



A(k, -k-p,p)A(k)D(k + p). 



(3.15) 



Once the above connection has been established, we return to the derivation of the explicit 
expression for the form factor A in the soft ghost limit. To that end, we will consider again 
the diagrams shown in panel (B) of Fig. HJ with dressed gluon and ghost propagators, and 
tree-level values for all the interaction vertices. In this configuration, the expressions given 
in Eq. (13. 2p reduce to 

(<*!)„ = p p {k + p)° J(l+p) v D(l+p)D(l + k+p)A(l)P po (l), 

(<fc)„ = P p {k+pf jD(l+p)A(l)A(l-k)P^l)P^l-k)T^ a , (3.19) 

We next outline the general procedure for isolating the A(— k, 0, k) defined in Eq. (13.131) . 
First, we observe that the most general Lorentz decomposition of the diagrams given in 
Eq. (I3TT91 is 

{di) v = p p (k + p) a [hg vp K + figuakp + f s g pa k v + hg vp p<j + hg^Pp + Ug P aPu 
+ fiVvPpPa + hPuPpK + fgPvkpPa + fioPukpK + fuKk p k a + f 12 k u k p p a 
+ fnKPpka + fuk v ppp a ] , (3.20) 

where the corresponding form factors /, = fi(k,p) are assumed to be finite in the infrared 
limit p — > 0. 

A detailed look at this expansion reveals that only the tensorial structure g U pk a , accom- 
panying the form factor fi, can saturate the pref actor p p {k + p) a and survive when the limit 
p — > is taken. Specifically, we may rewrite Eq. (13.201) as follows 

{di) v = p p ^f l {k,p)g up K + 0{p){k + p) u 

= k 2 f 1 (k }P )p u + 0(p)(k + p) U} (3.21) 

where the symbol 0{p){k + p) v is used to indicate terms that saturate with p u or k u , but 
whose form factors are of order 0(p) or higher, and will not contribute in the soft ghost 
configuration. Furthermore, one can perform the Taylor expansion of fi{k,p) around p = 0, 
namely, 

h(k,p) = fx(k,0) + 2(k ■ p)f[(k) + 0(p 2 ); = A/^p)]^. (3.22) 

Thus, only the zero order term of this expansion is relevant for our kinematic configuration, 
and we obtain finally from Eq. ( 13 . 2 1 [) the following result 

{di) v = k 2 h(k, OK + 0(p)(k + p) v , (3.23) 
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where the quantity k 2 fi(k,0) should be identified as the contribution of the corresponding 
diagram to A(k), while terms containing the derivatives of f\ are naturally reassigned to 
G(p)(k+p) u . 

After these observations, it is relatively easy to establish that this generic procedure can 
be systematically implemented by performing the following steps: (i) Set p = from the 
beginning inside the integrals of Eq. (I3.19p . (ii) Discard all the terms that give rise to 
structures of the type 0(p)(k + p) v . (iii) Determine the contribution of the diagram that 
saturates the index of the momentum p p with the metric tensor g vp . 

To illustrate in some detail the above procedure, let us focus our attention on the con- 
tribution of diagram (d\), appearing in the first line of Eq. ( 13.1 9p . Applying step (i), we 
obtain 

{d x ) v = p p (k + pY[ Jl v D{l)D{l + k)A(l)P pa (l) + p u Jd{1)D{1 + k)A(l)P pa (l) } . (3.24) 

Now, using criterion (ii), it is easy to recognize that the part of Eq. (13. 24ft to be retained is 
given by 

{d x ) v = -p p l vp {k) , (3.25) 

where we have defined the integral 

h P (k) = J^jP-D{l)D(l + k)A(l)l v l p , (3.26) 
which may be further decomposed as 

I vp {k) = h(k 2 )g up + h{k 2 )k v k p , (3.27) 

with 

Jl ^ 2 ) = dTI^WpW; J 2(A; 2 ) = ^_ (dW - k 2 g^)I up (k) . (3.28) 

Thus, using Eqs. (13.27P and (13.28p . we obtain from Eq. (I3.25j) the following result 

{di)u = -^Z\Vu \ { ^T^ - (l ■ k) 2 )D(l)D(l + k)A(l) , (3.29) 

where, according to (iii), we have only written explicitly the contribution that saturates the 
momentum p p with the metric tensor g vp . 
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Consider finally the contribution of diagram (d 2 ). After the shift I i— > —I, and setting 
p = inside the integral, it becomes 

(d 2 ) v = p'\k + pf Jd(1)A(1)A(1 + k)P;(l)P«(l + k)Yf aa . (3.30) 
It is then elementary to show that 

p»(k+pfP°(l)P^l + A;)rjl = 2 P P / 2(/ %2 ^ fc2 " ( l ■ k ? + V + k ) 2 ( l ■ k )\ 

+ 2p J l -^ + -l 2 /K o(p)(k + p) u: (3.31) 
and, therefore, the part of diagram (d 2 ) to be saved is 

{d 2 ) u = 2p u J [K > + l D(l)A(l)A(l + k)+ 2fQ up {k) , (3.32) 
where we have defined the integral 

Qvp{k) = / / 2 J+ P fc)2 [ /2fc2 - V ■ k ) 2 + + k ^ 1 ■ k)]D(l)A(l)A(l + k) . (3.33) 

One observes at this point that the first term in Eq. (13.321) is already saturated by p v 
and may be assigned to the form factor A(k) without further considerations. On the other 
hand, decomposing the integral Eq. (I3.33f) in the second term as 

Qu P (k) = Q x {k 2 )g vp + Q 2 (k 2 )k u k p , (3.34) 

with 

Qi(k 2 ) = ^-jP» p (k)Q up (ky, Q 2 (k 2 ) = — L— (d^ - k 2 g^)Q up (k) . (3.35) 

we obtain from Eq. (I3.32p the result 

(da)* = P ^fe 2 "^ ^f. fe^g 2 - 1 K^H-A;) 2 C^-fc) — C^-fc) 2 — C^— 2)^ 2 fc 2 ]^C0^C0AC/-»-A;) , (3.36) 

where, as before, we have omitted terms of the type 0(p)(k + p) v - 

Once Eqs. (I3.29|) and (13.361) have been derived, we will use Eq. (13.61) for projecting out 
the form factor A(k), as well as Eqs. (13.101) and (13.111) . in order to pass to Euclidean space, 
and subsequently cast the answer in spherical coordinates. Thus, we arrive at the final result 

(3.37) 



A(y) = / dty/ylF{t)A{t) d6' sin 4 9' cos 9' 

127T J Jq 
/~i /*oo pir 

+ < f±A / dtF(t)A(t) / cifl'sin 4 ^ 
6tt j Jo Jo 



u 

[yt(l + sin 2 9') - (y + t)^yi cos 9'} 



u 

Notice that, in this case, y = k 2 , u = (I + k) 2 , and 9' is the angle between k and I 
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Ghost dressing function 

a Lattice V=64 4 and (J=5.7 
4> Lattice V=80 4 and p=5.7 
Fit 




q [GeV ] 



FIG. 5: Lattice results for t he g luon propagator, A(g), (left panel) and ghost dressing, F(q), 



(right panel) obtained in Ref. 
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and renormalized at \x = 4.3 GeV. The (red) continuous curves 



represent the corresponding fits for the lattice data. 



IV. NUMERICAL RESULTS 



In this section we will carry out a detailed numerical analysis of the equations obtained 
in the previous sections. Specifically, in the first subsection we determine A(0,— p,p) by 
solving the integral equation Eq. (13. 12ft . using the lattice data of as input for the 
gluon propagator A(q) and the ghost dressing function F(q ) appearing in it. The solution 
obtained is then compared with the lattice data of 32j, [33J . In the second subsection, we 
solve numerically the coupled system formed by the integral equations of the ghost dressing 
function (12. 16[) and of the ghost-gluon vertex in the soft ghost configuration, given by (13.371) . 
The unique external ingredient used when solving this system are the lattice data for the 
gluon propagator A (q). The solution obtained for F(q) compares very favorably with the 
lattice data of 111 ]. 



A. Solution for the soft gluon configuration 



The integral equation (13.121) is solved through an iterative process, using as input for 
the gluon propagator and the ghost dressing function the data obtained from the SU(3) 
quenched simulations of [llj, shown in Fig. [51 Note that the lattice data shown have been 
renormalized at /J, — 4.3 GeV, within the MOM scheme. The value of a s that corresponds 
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to this value of \x may be obtained from the higher-order calculation presented in 38| ; 
specifically, we have that a a (fi) = 0.22. 



1.4 



1.3 



1.2 



Q. 



Z- 1.1 - 



1.0 



0.9 




Gluon-ghost form factor A^O.-p.p) 
Soft gluon configuration 

o Lattice V=32 4 and (S=5.8 
s Lattice V=32 4 and p=6.0 
* Lattice V=48 4 and p=6.0 
A^O.-p.p) for a(n)=0.22 



10 



p[GeV] 



FIG. 6: Numerical result for A(0, —p,p), obtained from Eq. (|3.12p when a s {(£) = 0.22. 



The (red) continuous line in Fig. [6] represents the corresponding solution for A(0, —p,p). 
We clearly see that A(0, — p,p) develops a sizable peak around the momentum region of 
830 MeV. In addition, as had been anticipated in the subsection IIIIB[ we confirm numeri- 
cally that A indeed assumes its tree level value when p — > 0, i.e., A = 1. It is also interesting 
to notice that, in the ultraviolet limit, the form factor gradually approaches its tree level 
value. 

In Fi g. El we compare our numerical results with the corresponding lattice data obtained 



in Ref. 



32 



for this particular kinematic configuration. Although, the error bars are 
rather sizable, we clearly see that our solution follows the general structure of the data. In 
particular, notice that both peaks occur in the same intermediate region of momenta. Evi- 
dently, A(0, —p,p) receives a significant non-perturbative correction, deviating considerably 
from its tree level value. 
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FIG. 7: Left panel: The form factor A{— k, 0, k) (circles) and the fit given by Eq. (|4.3p (red 
continuous line). Right panel: The numerical solution of F(p) (red continuous line) compared with 



the lattice data of Ref. 



111 ]. Note that the value of a s used when solving the system is a s (fi) = 0.22. 



B. The coupled system: ghost SDE and ghost-gluon vertex. 

In this subsection we present the central result of the present article, namely the modifi- 
cations induced to the ghost dressing function by the inclusion of a non-trivial structure for 
the corresponding ghost-gluon vertex. 

To that end, after passing to the Euclidean space and introducing spherical coordinates, 
using Eq. (I3.10P and Eq. (13. lip , we obtain from Eqs. (I2.16P and (I3.37P the expressions 



F-\x) 



2tt 2 



dyyA(y)A(y) / d6sm A 6 
(I Jo 



F(z) F(z') 



(4.1) 



and 



Mv) 



i - 



a s C A 



+ 



127T 2 

a s C A '" x 



POO pTT 

I dty/ylF(t)A(t) / d0'sin 4 0'cos0' 
o Jo 



F(u) 



u 



(4.2) 



6tt 2 



O P7V 

dtF(t)A(t) / c^'sinV 
Jo 



A(u) 



u 



[yt(l + sin 2 0') - (y + t)y/yi cos 9'} 



where now z = (k + p) 2 , z' = (k + fi) 2 and \x is the renormalization point introduced within 
the MOM scheme, i.e., by requiring that F' 1 ^ 2 ) = 1. 

We next solve the above system iteratively, using again the lattice data for A(q) and 
a s (n) = 0.22 as input. The results for F(p) and A(—k, 0, k) are shown in Fig. [71 

On the left panel of Fig. [7J the curve in circles represents the result for A(—k,0,k). 
Evidently, A develops a peak in the intermediate region of momenta, in a way similar to 
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the case discussed in the previous subsection. In this case the maximum of the peak occurs 
around 1 GeV, and once more, in the infrared and ultraviolet limits A(—k, 0, k) assumes its 
tree-level value. 

On the same panel we show a fit for A{— k, 0, k), represented by the (red) continuous 
curve, whose functional form is given by 



ak 2 

[(k 2 + b) 2 + c\\n (d + k 2 /k 2 ) 



A{-k, 0, k) = 1 + r/ ,„ 2 | u , 2 | | ,„ 2/ ,„ 2 , , (4-3) 



with the following values for the fitting parameters a = 0.68 GeV 2 , b = 0.72 GeV 2 , 
c = 0.29 GeV 4 , d = 9.62 and k 2 = 1 GeV 2 . 

On the right panel of Fig. 0, we compare our numerical result for F(p) (red continuous 
curve) with the corresponding lattice data of Ref. llj, observing a rather notable agreement. 
We emphasize that, contrary to what happens when the bare vertex is used (see Fig. [3]), 
the accuracy achieved here does not rely on the artificial enhancement of the value of the 
coupling; the latter, as mentioned above, was kept at its standard value predicted from 
general MOM considerations. 

It is important to realize that, although A does not provide a sizable support for ghost 
SDE in the deep infrared, the contribution that it furnishes in the region of intermediate 
momenta is sufficient for increasing the saturation point from F(0) = 1.67 to F(0) = 2.95 
(Figs. [3]and[7J respectively). This observation suggests that the ghost SDE is particularly 
sensitive to the values of its ingredients at momenta around two to three times the QCD 
mass scale. 



V. CONCLUSIONS 



In the present work we have considered the "one-loop dressed" approximation of the 
SDE that governs the evolution of the ghost-gluon vertex. In particular, we have focused on 
the dynamics of the form factor denoted by A, which is the one that survives in the SDE 
for ghost dressing function, in the LG. The vertex SDE has been evaluated for two special 
kinematic configurations, one of them corresponding to the well-known Taylor limit. When 
coupled to the SDE of the ghost, the contribution of this particular form factor accounts for 
the missing strength of the associated kernel, allowing one to reproduce the lattice results 
rather accurately, using the standard value of the gauge coupling constant. 
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The fact that, despite the truncation implemented on the vertex SDE, we finally obtained 
a rather good agreement with the lattice, hints to the possibility that the omitted terms are 
numerically subleading, at least in the case of the special kinematic configurations consid- 
ered. It might be interesting to pursue this point further. Specifically, in the present analysis 
the terms proportional to the second form-factor, denoted by B, have been automatically 
discarded, precisely because they do not contribute to the ghost SDE. However, given that 
both form factors participate in the fundamental relation of Eq. (|2.6p . one might consider the 
possibility of keeping these terms throughout the calculation, and then checking explicitly 
to what extent Eq. (I2.6P is satisfied in the present approximation. 

Recently, the study of the effects that the dynamical quarks induce on some of the 
fundamental Green's functions of QCD has received particular attention, both from the 
point of view of unquenched lattice simulations Q|, as well as by means of an SDE-based 
approach 41]. In particular, lattice simulations reveal that the inclusion of light active 
quarks results in a considerable suppression in the deep infrared and intermediate momentum 
region of the gluon propagator. This characteristic feature has been firmly established also 
within the SDE framework of 4l|. On the other hand, the unquenched ghost dressing 



function simulated on the lattice suffers minimal changes from the inclusion of quarks [401 ] : 
this property has also been anticipated within the aforementioned SDE analysis 41], as a 
direct consequence of the fact that, in the case of F, the quark- loops enter as "higher-order" 
effects. In addition, it is well-known that the value of the MOM coupling, a(/x), increases 
in the presence of quark loops. 

It would be, therefore, interesting, to study the combination of these competing effects 
systematically, including the vertex equation for A, derived here. In particular, the non- 
linear nature of the corresponding integral equations converts this combined analysis into 
a rather challenging problem. Specifically, the changes induced to the integral equation for 
A, due to the aforementioned suppression of the gluon propagators entering in it, must be 
compensated, to a considerable level of accuracy, by the corresponding increase in the cou- 
pling constant, in order to finally obtain the rather minor change observed in F. We hope 
to be able to carry out such a study in the near future. 



21 



Acknowledgments 



The research of J. P. is supported by the Spanish MEYC under grant FPA2011-23596. 
The work of A. C. A is supported by the National Council for Scientific and Technological 
Development - CNPq under the grant 306537/2012-5 and project 473260/2012-3, and by 
Sao Paulo Research Foundation - FAPESP through the project 2012/15643-1. 



[1] C. D. Roberts and A. G. Williams, Prog. Part. Nucl. Phys. 33, 477 (1994). 

[2] R. Alkofer, L. von Smekal, Phys. Rept. 353, 281 (2001). 

[3] C. S. Fischer, J. Phys. G G32, R253-R291 (2006). 

[4] D. Binosi and J. Papavassiliou, Phys. Rept. 479, 1-152 (2009). 

[5] A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D 65, 025012 (2002). 

[6] A. P. Szczepaniak, Phys. Rev. D 69, 074031 (2004). 

[7] A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. D 78, 025010 (2008). 
[8] A. C. Aguilar and J. Papavassiliou, Phys. Rev. D 83, 014013 (2011). 

[9] A. Cucchieri and T. Mendes, PoS LAT2007, 297 (2007); Phys. Rev. Lett. 100, 241601 (2008); 

Phys. Rev. D 81, 016005 (2010); PoS LATTICE2010, 280 (2010). 
[10] P. O. Bowman, U. M. Heller, D. B. Leinweber, M. B. Parappilly, A. Sternbeck, L. von Smekal, 

A. G. Williams and J. -b. Zhang, Phys. Rev. D 76, 094505 (2007). 
[11] I. L. Bogolubsky, E. M. Ilgenfritz, M. Muller-Preussker and A. Sternbeck, PoS LATTICE, 290 

(2007). 

[12] I. L. Bogolubsky, E. M. Ilgenfritz, M. Muller-Preussker and A. Sternbeck, Phys. Lett. B 676, 
69 (2009). 

[13] O. Oliveira, P. J. Silva, Phys. Rev. D79, 031501 (2009). 

[14] O. Oliveira and P. J. Silva, PoS LAT2009, 226 (2009). 

[15] M. R. Pennington and D. J. Wilson, Phys. Rev. D 84, 119901 (2011). 

[16] A. Bashir, L. Chang, I. C. Cloet, B. El-Bennich, Y. -X. Liu, C. D. Roberts and P. C. Tandy, 

Commun. Theor. Phys. 58, 79 (2012). 
[17] P. Boucaud, J-P. Leroy, A. L. Yaouanc, J. Micheli, O. Pene and J. Rodriguez-Quintero, JHEP 

0806, 012 (2008). 



22 



[18] T. Kugo and I. Ojima, Prog. Theor. Phys. Suppl. 66, 1 (1979). 

[19] T. Kugo, |arXiv:hep-th/9511033| 

[20] P. Watson and R. Alkofer, Phys. Rev. Lett. 86 (2001) 5239. 

[21] K. -I. Kondo, Phys. Rev. D 84, 061702 (2011). 

[22] D. Dudal, J. A. Gracey, S. P. Sorella, N. Vandersickel and H. Verschelde, Phys. Rev. D 78, 

065047 (2008). 

[23] J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 

[24] A. C. Aguilar and J. Papavassiliou, JHEP 0612, 012 (2006). 

[25] A. C. Aguilar and A. A. Natale, JHEP 0408, 057 (2004). 

[26] D. Binosi and J. Papavassiliou, Phys. Rev. D 77(R), 061702 (2008). 

[27] A. C. Aguilar, D. Binosi and J. Papavassiliou, JHEP 1201, 050 (2012). 

[28] D. Dudal, O. Oliveira and J. Rodriguez-Quintero, Phys. Rev. D 86, 105005 (2012). 

[29] A. C. Aguilar, D. Binosi and J. Papavassiliou, JHEP 1007, 002 (2010). 

[30] A. Cucchieri, T. Mendes and A. Mihara, JHEP 0412, 012 (2004). 



[31] E. M. Ilgenfritz, M. Muller-Preussker, A. Sternbeck and A. Schiller, [arXiv:hep-la t/0601027 
[32] E. -M. Ilgenfritz, M. Muller-Preussker, A. Sternbeck, A. Schiller and I. L. Bogolubsky, Braz. 
J. Phys. 37, 193 (2007). 



[33] A. Sternbeck, hep-lat/0609016 



[34] A. Cucchieri, A. Maas and T. Mendes, Phys. Rev. D 77, 094510 (2008). 

[35] W. Schleifenbaum, A. Maas, J. Wambach and R. Alkofer, Phys. Rev. D 72, 014017 (2005). 

[36] J. C. Taylor, Nucl. Phys. B 33, 436 (1971). 

[37] W. J. Marciano and H. Pagels, Phys. Rept. 36, 137 (1978). 

[38] P. Boucaud, F. De Soto, J. P. Leroy, A. Le Yaouanc, J. Micheli, O. Pene and J. Rodriguez- 
Quintero, Phys. Rev. D 79, 014508 (2009). 

[39] A. C. Aguilar, D. Binosi, J. Papavassiliou and J. Rodriguez-Quintero, Phys. Rev. D 80, 085018 
(2009). 

[40] A. Ayala, A. Bashir, D. Binosi, M. Cristoforetti and J. Rodriguez-Quintero, Phys. Rev. D 86, 
074512 (2012). 

[41] A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. D 86, 014032 (2012). 



23 



